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We predict the rates of the charmless three-body 5° — * p~pn and TT~pn modes due to weak vector 
current contributions to be 4 x 10~^ and 2 x 10~®, respectively. The basis is a factorization 
approach of current produced nucleon pairs, together with an isospin transformation that relates 
nucleon weak vector form factors to electromagnetic form factors. Adding the axial vector current 
contribution, we find p~pn and B'^ — > p^pn to be at 10~^ order. The three-body modes 

appear to dominate over the two-body modes such as _B — > pp, pA. 
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PACS numbers: 13.25.Hw, 13.40. Gp, 14.20.Dh 
I. INTRODUCTION 

A large number of charmless mesonic decays of the B 
mesons have emerged since 1997, and are of great current 
interest. For example, vr+vr^/A'+Tr^ ~ 1/4— 1/5 sug- 
gests [0 that <^3 (or 7) = arg V*^^ could be 90° or more. A 
natural question to ask Q is: what about charmless bary- 
onic decays? The CLEO Collaboration has done some 
search in the past, but turning up null results for modes 
like B° pp that are below the 10^^ level [|. The 
Belle Collaboration has recently improved the limits ^] 
on B'^ — > pp by an order of magnitude, pushing down to 
the 10^^ level. It is of interest to ask, therefore, if all 
charmless baryonic modes are below 10~^. In this paper, 
from a suitably sound theoretical basis involving nucleon 
form factor data, we show that _B" p~pn could be a 
leading charmless baryonic decay with rate at the 10~^ 
level. 

The CLEO Collaboration recently reported the obser- 
vation of the B'^ — > D*~pn mode at the 10~^ level [||, 
which is only a factor of 4-5 lower than B° D*^ and 
I?*~7r+ Scaling by |T46/KbP one could already infer 
that B'^ p~pn ~ 10~^, but a better understanding 
is desirable. A factorization approach for B^ D*^pn 
with current produced nucleon pairs has been proposed 
recently [Q. The three-body decay is seen as generated 
by two weak currents: one converting to D*~ , the 
other creating the nucleon pair. The nucleon weak vec- 
tor form factors are related by isospin rotation to nucleon 
electromagnetic form factors. By using G^" measured 
from e+e" NN and pp e^e~ processes p[-[Il], we 
are able to account for up to 60% of the observed rate, 
the remainder seemingly coming from axial vector cur- 
rent contribution. Emboldened by this success, we apply 
the approach to the charmless modes B'^ p pn, ■n^pn 
where one replaces the by p or tt. 

Besides rates, we are able to predict the pn pair mass 
spectrum. Since the vector current contributes domi- 
nantly to the total rate for the D*~pn case Q, we 

expect the vector current to dominate the — > p~pn 
and TT~pn rates as well. Incorporating estimates of the 
axial current contributions, the total rates are slightly 
higher than from the vector current alone. 



II. FORMALISM 

Our starting point is to factorize the current produc- 
tion of ph pairs, i.e. 



{p-{TT-)pn\neff\B") = :^V,dV:,a, 



Gp 

7f 



x{p-{^-)\V^ - A^\B''){pn\V^ - A^\ 0) . (1) 

The V — A current induces the b u transition in the 
first term, while in the second it creates the nucleon pair, 
which is illustrated in Fig. |^. It is an extension of the 
usual factorization of current-current matrix elements in 
the case of B decays to two mesons [|2|,^ . 

Once factorized, the matrix element {pn\Vf^ — ^;^|0) 
describes nucleon pair creation by a charged weak cur- 
rent. The matrix element for the vector (V^) portion 
can be expressed as 



{pn\V+\0) ^ uip,) 



Frith, +^ 



2mN 



(2) 



where m^r is the nucleon mass, q = {pp + Pn) the mo- 
mentum transfer, t = q^ = rrip^ the p-n invariant mass 
squared, and i^j^ are nucleon weak form factors, with 
F^it) normalized at t = Q, 



(0) = 1. 



(3) 



The photon field contains , which, together with 
form a weak isotriplet. The coupled currents also 
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FIG. 1. Feynman diagram illustrating Eq. 
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form an isotriplet, and can be interrelated via an isospin 
transformation. For the nucleon pair, the strong isospin 
symmetry of the nucleon state coincides with the weak 
isospin symmetry of the weak and em currents. The weak 
vector form factors are therefore related to electromag- 
netic (em) isovector form factors. 

The matrix element {N(j)')N{p)\Jli'^\Q) for the em 
current can be expressed as 

{N{p')N{p)\j;nO) 

= u{p')^^F^{t)^, + r ^a,.q^y{p), (4) 

where i^i,2(i) are respectively the Dirac and Pauli form 
factors, normalized at t = as 

i^f(0) = l, i^r(0)=0, i^fW^Kp, F2"(0) = K„, 

(5) 

with ftp (•„■) the proton (neutron) anomalous magnetic 
moment in nuclear magneton units. These form factors 
are related to the Sachs form factors via 



t 



GTit) = Frit) + 



Frit): 



Grit)^Frit)+Frit). (6) 

The isospin decomposition of the em current is given by 
1 



{F[±Fn.. 



i = 1,2, 



(7) 



where s, v stand for the isoscalar and isovector compo- 
nents, respectively. The isovector component of the em 
current and the vector portion of the charged weak cur- 
rents form an isotriplet, as manifested by |14| 

2F,^t)^Frit), 



z = l,2. 



(8) 



where the factor 2 is from the definition of Fi'^2^\t) in 
Eq. For example, from Eqs. (^), (||) and (0) one 

easily checks that 2F|'(0) = F^{0). 

We can now write the three-body — > p~pn decay 
amplitude in the following form 



■ GP TA* 



-ig^^u{mB +mp)Ai{q^) + i{pB +Pp) 



xu{pp) 
and for i?° 



2iF^+F^)r + ^{p,,^p,r 

TUN 



rriB + TTLp 
A2{q^) 

v{Pn), (9) 



TT pn, 



G 



iMp = -i —^VudV*f, ai {ps +Ptt),, Fi{q 



V2 



xu{pp) 



F- 



2{F^ + F^)r + — {P^ 
rriN 



Ppf 



v{Pn), (10) 



TABLE I. Form factors at q'^ — and the parameters ai.; 
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0.73 


1.40 


0.48 
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0.50 





where Cp is the p meson polarization, and V{q'^), Ai{q^), 
^2(9^) and Fx{q-^) (not to be confused with baryon form 
factors) are the transition from factors arising from the 
6 — !■ u transition in the first matrix element of Eq. (|l|) . 

As we need to integrate over for these three-body 
decay modes, we need to pay more attention to the q^ 
dependence of these transition form factors ||l5|. How- 
ever, since our focus is on utilizing experimental data on 
baryon form factors, for _B — > p, tt form factors, we shall 
simply take what is readily available in the litetature. 
Among the several recent models for the meson form fac- 
tors (see e.g. Ref. |p^ , |l7| ), we shall use |17 



/(O) 



{l-q^/Ml){l-a^q^|M^y 
for i^i(g2) and V{q^), and 

/(O) 



lie) = 



l-a^q-^/M^ + a^qVMl 



(11) 



(12) 



for Ai 2(9^). My is the appropriate pole mass which is 
taken to be 5.32 GeV. Note that the q^ dependence is 
quite different from the monopole form used in Ref. . 
In Table || we give the values of the relevant form factors 
at zero momentum transfer as well as the parameters 

It is important to note that the baryon form factors 
must satisfy perturbative QCD (PQCD) quark counting 
rules |l^, which give the leading power large- |i| fall-off 
of the Fi{t) form factor. Since helicity-flip gives an extra 
1/t factor for F2{t), one finds in the large \t\ limit 



i-1,2, (13) 



where Qo ~ Aqcd = 0.3 GeV, 7 = 2-1- 4/(3/3), and 
(3 is the QCD /3-function to one loop. We note that 7 
depends weakly on the number of flavors; for three fia- 
vors 7 = 2.148. The asymptotic form given in Eq. ( |T^ ) 
has been confirmed by many measurements of the nu- 
cleon form factors = Ff'" + Fr over a wide range 
of momentum transfers in the space-like region [p^ . It 
has also been confirmed in the time-like region with the 
recent nucleon time-like data ||9|-[ll|] . 

The combination 2{F^ -f ^2") in Eq. ^ can be re- 
placed by G\j — G^ which is composed of measurable 
quantities. Similar replacement can also be made for 
F2, which is a combination of — and G^^ — G^. 
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Most time-like data for the magnetic form factors, how- 
ever, are extracted by assuming either \G^\ = \G^j\ or 
in the explored region of momentum transfer. 

= (1 - t/Am%) clearly vanishes at 
threshold, by assuming \G^\ — \G^i\ in extracting GjJ^ 
from data, the information on is lost. In our calcula- 
tion we concentrate on the part of Eq. (S) which contains 



Since G^j 



pv 



F2 , the contribution from F2 can be determined 



only when Gj^ and G^ can be separated from data with 
better angular resolution. 

We take \G^.j\ in the following form [Q to make a phe- 
nomcnological fit of the experimental data |^-|Tl|] : 



GUt)\ 



5 

i=l 



In 



(14) 
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FIG. 2. Time-like proton magnetic form factor data, fitted 
by Eq. (jl^) with the parameters given in Eq. (|l6|). 



(15) 



the power of the leading term and logarithmic factor are 
as suggested by PQCD, and the fewer number of fit pa- 
rameters for G^^ reflects the fact of scarcer neutron data. 
We find the best fit values 

xi^ 429.88 GeV", - - /f/focoo r.„^rio 

-10783.69 GeV^ 



= -448583.96 GeV 
X5 = 635695.29 GeV^^ 



X3^ 109738.41 GeV 
and 
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(16) 
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□Antonelli et al 9B(|GeI=0) 
■ Antonelli et al 98 ( | Ge I = I Gm I ) 
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yi = 236.69 GeV* 



y2 



^579.51 GeV'^ 



(17) 



where the per degree of freedom (d.o.f.) of the fits are 
1.39 for |G^,^| and 0.41 for \Gli\, respectively. We show 
in Figs. ^ and |^ the best fit curves given by Eqs. ([l^) and 
( p^ ) with the above parameters. 

It was pointed out in Ref. [|o| that the data supports 
|G^| = as well. We therefore perform a fit to the 
neutron magnetic form factor data extracted under the 
0, giving the best fit values 



assumption of |G^ 



yi = 292.62 GeV'' 



2/2 = -735.73 GeV^ 



(18) 



with x^/d.o.f.= 0.39, which is slightly lower than the 
previous fit. This fit is also plotted in Fig. ^. More data 
is needed to distinguish between the two cases. 

We note that there is a sign difference between G^^ 
and G^j in the space-like region. Since analyticity implies 
continuity at infinity between space- like and time-like pc| ] 
regions, the time- like magnetic form factors are expected 
to have similar behavior as the space-like ones: real and 
positive for the proton, but negative for the neutron. 

For large t, QCD predicts the magnetic form factors 
to be real [p^ , with the neutron form factor weaker than 
the proton case ||2l|. According to QCD sum rules [ p2[ , 
asymptotically one expects G\j/G\,j ^ Qd/Qu = —0.5. 



FIG. 3. Time-like neutron magnetic form factor, where 
the solid (long-dash) line is the fit given by Eq. ( |l5| ) with 
parameters given in Eq. (|l^) (Eq. (^3)), for data extracted 
with \Gl\ ^ \Gli\ {\G"e\ = Q) assumption. 



In our fits, with sign difference between G\j and G"/, 



we have GlJGlj = -yi/xx = -0.55 (-0.68) for \Ge\ 
\Gm\ {\G%\ — 0). Nucleon form factors have also been 
analyzed from negative to positive t with dispersion re- 
lations. The phase of the proton magnetic form factor 
turns out to be 27r, hence the proton magnetic form 
factor is real and positive as expected asymptotically, 
starting already from t > 4 GeV^ P^P^ onwards. 



III. RESULTS AND DISCUSSION 

We still need to fix VudV*^ai. We shaU take = 
0.9747 and |Kb| = 34.95 x lO"'^ from Ref. ||. For the 
effective coefficient ai, we take the value ai = 1.05 from 
Ref. ||l^ for effective number of color = The sit- 
uation is slighly different from previous study on — )■ 
D*^pn case |^], where oi is taken from D*~ de- 

cay ||l^ . Here, when one tries to follow the procedure by 
looking at — > P~ 7r~p+ modes, the tree-penguin 
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FIG. 4. The vector current induced differential de- 
cay rates dVv{B° ^ p~pn)/dq^ (upper two curves) and 
dVv{B° — > 7r~prl)/dg^ (lower two curves). Solid lines are 
from fitting nucleon form factor data with |G5/"I ~ 
dashed lines with IG^J = IG^^I and |G^.| = 0. 



fast recoil meson carries away energy and would be more 
favorable for baryon production in the recoil system be- 
cause of reduced energy release compared to two-body 
decays. This fast recoil meson accompanying the low 
mass baryon pair can be tested experimentally. We note 
that the narrowness of the t distribution in Fig. ^ implies 
the recoil meson spectrum of a quasi-two-body mode. 

According to the i?° D*~pn study [0, the vector 
current contributes ~ 60% of the observed rate |^ . We 
can estimate the total rates of the P~pn, TT~pn 

modes by assuming similar proportions of the vector 
current contributions. The estimated rates are then 
Bt{B° p-pn) - 7 x 10^6 and Br(7r-pn) ~ 3 x 10-^ 
i.e. B'^ p^pn is of order 10^^. Inspection of Fig. 1 
suggests that — > p^pn is at the same order. 

It is interesting to compare the three-body rates 
Br(_B" hpn) with the two-body ones Br(_B" hp'^) 
where h stands for the recoil meson. We find a similarity 
of the ratios Br(3-body)/Br(2-body) between h — D*^ 
and h = p^.n^. By taking 03 — 54.8° |25|, we obtain 
Bi{p-p+{TT'p+)) ~'32 (22) x 10-6 We then have 

interference will complicate things 0,|2^. We therefore Bi{p~pn{'K~pn))/Bi{p~ p'^{n^ p^)) ^ 0.22 (0.14) which 

is rather close to BT{D*-pn)/BT{D*- p+) ~ 0.2. 

So far, we have assumed only the tree level b — > duu 
transitions as the underlying process, as illustrated in 
Eq. (1) and Fig. 1. As we make comparison with charm- 
less mesonic modes, it is important to remember that, 
unlike the B^ — )■ D*~pn case, we expect significant pen- 
guin contributions as well. It is known that b — > suu 
penguins dominate Kt: modes, and that penguin contri- 
butions make significant impact on tttt modes. Since the 
impact of penguins becomes less pronounced for VP and 
VV modes, and since our approach to p~pn corresponds 
to both VP and VV components, we expect the impact 
of penguins on — )■ p^pn to be not as pronounced 
as in i?*^ — !■ vr+TT- . On the other hand, extending our 
discussion to 6 — > s penguin operators suggests three- 
body charmless baryonic modes such as p~Ap, K*A'S 
etc. However, this would involve further assumptions. 

With recent improvement of i?*^ — > pp limit to 10"^ 
level ji] , we now have the intriguing situation that B — > 
ppn ^ B ^ pp, i.e. the three-body baryonic mode dom- 
inates over the two-body. This would confirm the con- 
jecture made in Ref. We urge the Belle and BaBar 
groups to search for the three-body modes experimen- 
tally. 

In summary, from a relatively robust foundation, in 
analogy with the _B" D*~pn mode observed by CLEO, 
we have shown that i?° — > p^pn is very likely at the 10"^ 
level and should be well within the capabilities at B Fac- 
tories. Together with the absence of two body modes 
such as B^ —> pp, the dynamics of B decay to baryons 
are somewhat different from mesonic final states. Ex- 
tending our study to include the standard set of effective 
operators should be straightforward, and one expects a 
host of three-body baryonic modes at the 10"^ level. We 
expect charmless baryonic modes to emerge soon at the 
B Factories, likely from three-body onwards. 



use the short distance ai for simplicity. 

For proton and neutron data extracted assuming 
— I^m"Ij find the branching ratio for B^ 
p~pn arising from the vector current is 

BMB° ^ p-pn) = (3.58tS]?) x lO'^ (^)' , (19) 

where the subscript y is a reminder that this is from 
the vector portion of the weak current alone. The upper 
and lower bounds correspond respectively to the maxi- 
mum and minimum of the branching fraction evaluated 
by scanning through ^ Xmin + 1 in the fits. In a sim- 
ilar fashion, for data extracted assuming |G^| = |G^/| 
but \G'^\= 0, we find 

Brv(B" ^ p-pn) = (4.53tJ:°^) x lO'^ (^)' . (20) 

Following the same methods, we also give predictions 
on the B^ — s- ir-ph rate that arise from the vector cur- 
rent. We find 

Btv{B' ^ n-pn) = (1.82^^11) x (^)' , (21) 



for \GP^\ 



IG^I, and 



Brv(i3° - n-pn) ^ (2.29l°:l^) x IQ-' (^)' , (22) 

for \Gyi = \GU hut |G1;|=0. 

Fig.^ shows the vector current induced differential de- 
cay rates dTv/dq^ of both the p-pn and the Ti-pn 
modes with BSW form factors. The peaking of the dif- 
ferential rates at ~ 5 GeV^ and ~ 4 GeV^ for the p-pn 
and n-pn modes, respectively, is a threshold enhance- 
ment effect for baryon production. As argued in 0, a 
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